Formation of an Edge Striped Phase in Fractional Quantum Hall Systems 
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We have performed an exact diagonalization study of up to N = 12 interacting electrons on a disk 
at filling v = | for both Coulomb and V\ short-range interaction for which Laughlin wave function 
is the exact solution. For Coulomb interaction and N > 10 we find persistent radial oscillations in 
electron density, which are not captured by the Laughlin wave function. Our results srongly suggest 
formation of a chiral edge striped phase in quantum Hall systems. The amplitude of the charge 
density oscillations decays slowly, perhaps as a square root of the distance from the edge; thus the 
spectrum of edge excitations is likely to be affected. 

PACS numbers: 73.40.H, 71.10.Pm 
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An important open problem in the physics of the frac- 
tional quantum Hall effect (FQHE) jjJD is the structure 
and the properties of edge states. The bulk FQHE is 
reasonably well understood: the kinetic energy of 2D 
electrons is quenched by the strong perpendicular mag- 
netic field B, and the Coulomb interaction dominates 
the physics of the partially filled Landau level (LL). At 
certain filling factors v the electrons condense into highly 
correlated gapped quantum fluids, which results in quan- 
tization of the Hall conductance a xy and vanishing diag- 
onal conductivity a xx There has been much interest 
in the physics of FQH edge states for several reasons. 
First, on a QH plateau the bulk electron states are local- 
ized by disorder, and gapless excitations are possible only 
in extended edge states ||||; thus the transport current 
flows along the edges . Second, FQH edge states have 
played an important role in many experiments: inter- 
nal resonant tunneling from one edge to another through 
states bound on a quantum antidot has been employed 
to measure the fractionally quantized charge of the tun- 
neling quasiparticles [0, and external electron tunneling 
into an edge of a FQH system has been used to study 
its excitation spectrum ||. Also, an edge channel is 
usually assumed to be localized within a few magnetic 
lengths £q = ^Jhc/eB at the boundary of the FQH sys- 
tem; then the low-energy dynamics is effectively ID, and 
field-theoretic descriptions of edge channels as a chiral 
Tomonaga-Luttinger liquid have been developed ||[l(J . 

Our quantitative understanding of the FQHE is based 
on Laughlin wave function 0] , which is known to be quite 
accurate for the bulk, translationally invariant v = | 
FQH state from extensive comparisons with exact di- 
agonalization results in the edgeless spherical geometry 
jfllO. However, the accuracy of this wave function re- 
mains relatively untested at the edge. Validity of a wave 
function in the bulk does not necessarily extend to the 
edge: the ground state of the system is separated from 
excited states by a gap in the bulk, which makes it in- 
sensitive to perturbations; on the other hand, the system 
has gapless excitations at the edge, which makes it more 
susceptible to the particular form of the interelectron in- 



teraction. The FQHE has been relatively less studied in 
the disk geometry |3, ^3| , H , because, in part, the inter- 
est was focused on the bulk states, and it is difficult to 
separate disk edge effects for small number of electrons. 

In this paper we report results of a detailed numeri- 
cal investigation of the microscopic structure of the FQH 
edge. To this end, we diagonalize the interaction Hamil- 
tonian in the disk geometry for up to N — 12 Coulomb- 
interacting electrons. For A > 10 we observe formation 
of a striped order in the charge density at the edges. 
These edge stripes are not captured by the Laughlin wave 
function that works well for the bulk FQH states. We also 
obtain an analytical fit to the numerical data, which sug- 
gests that the amplitude of the charge density oscillations 
decays only as a power law (specifically, as inverse square 
root) with distance, appreciably extending into the sam- 
ple. We interpret these results as formation of an edge 
striped phase (ESP) with wave vector g ESP w tt/2£q, pos- 
sibly smectic liquid crystal, at the edge of a FQH system. 

We study the simplest FQH state, that of spin- 
polarized electrons restricted to lowest Landau level, as 
appropriate in the large B limit, at filling v = |. The 
interaction Hamiltonian in the second-quantized form is 
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where cj„ creates an electron in the single-particle eigen- 
state of the angular momentum |m), where m = 0,1, .... 
There is no external confinement; the Coulomb matrix 
elements can be expressed via finite sums |l7| , |l8|] . 
The bulk ground state of N interacting electrons is well 
approximated by the Laughlin wave function: 
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N 
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where Zj = Xj — iyj (in units of £ ) gives the position of 
j-th electron as a complex number. This wave function 
is known to be the exact ground state |ll]|l5|] of the elec- 
trons interacting via short-range potential V% defined by 
expressing the Coulomb interaction as 



1 



V(\ Zj -z k \)=J2 V i P i jk . ( 3 ) 
i 

where Vi are the Haldane pseudo-potentials, and P^ k is 
the projection operator which selects the states in which 
particles j and k have relative angular momentum I [ fl6| . 

We construct \Pl and numerically as the ground 
states of the short-range and Coulomb Hamiltonians, 
respectively, for the total angular momentum M — 
^N(N — 1) for N electrons, which gives filling v = | 
in thermodynamic limit N — > oo. The Hilbert space 
is restricted by consideration of single-particle orbitals 
with angular momentum m < m max only. For V\ 
short-range interaction, the contribution of orbitals with 
m > TO^ ax = 3 (AT — 1) vanishes identically; for Coulomb 
interaction m^ ax is obtained by increasing Hilbert space 
untill overlap (^lI^c) converges to at least three signif- 
icant digits. For example, for N = 12, the largest v = | 
FQH system studied, M — 198 and the size of the Hilbert 
space is 15,293,119 for = 35. To the best of our 

knowledge, the largest v = | systems studied in the disk 
geometry prior to this study was N = 10 by Cappelli et 
al. |Q, who used only V\ interaction. 

In Fig. 1 we present the radial electron density profiles 
PL and pc for both 'J'l (short-range interaction) and 'E'c 
(true Coulomb interaction) for N = 5 to 12 electrons. 
A comparison of phif) with the exact Coulomb pc(?*) 
should clarify the importance of the long-range part of 
the Coulomb interaction for the FQH edge structure. For 
all AT, the exponential fall-off of pc{ r ) at the edge and 
the density oscillation nearest the edge is captured well 
by the Laughlin PlM- As expected, the edge shifts to 
larger r in accordance with r G d gc ~ (2mJ^ ax ) 1 / 2 . The 
properly normalized overlaps ^c) are given in Ta- 
ble I; the corresponding overlaps are much closer to unity 
in the spherical geometry, which indicates that the devia- 
tion originates in the disk edge. Still, given the huge size 
of the Hilbert space, the overlaps indicate that Laughlin 
'J'l capture the important short-range FQH correlations 
remarkably well, even at the edge. Table I also gives W, 
the weight of the components of *S?q with m > ni^ ax . 
All W < 1%, which means that fixing total angular mo- 
mentum M = ^N(N — 1) selects the v = ^ state for 
Coulomb interaction too: it fixes average density on the 
disk (pc) ~ (pl) (f° r r < ?~cdgc) to better than 10 -3 . 

For N = 6 to 8, density profile pc clearly displays a 
commensurability effect, strongest for N — 7, which can 
be visualized in a classical picture as tendency to have 
one central electron at r — surrounded by a ring of six 
electrons. For N — 7 the position of the outer maximum 
in the density profile is very close to the lattice constant 
a wc = (4tt % /3) 1/2 4 = 4.665£ , of the v = \ 2D Wigner 
crystal (WC) in the thermodynamic limit JL9|| . This com- 
mensurability effect is barely visible for the Laughlin den- 
sity profile pl (r) , indicating importance of the long-range 
part of Coulomb interaction for formation of an ordered 
phase, even for such small electron systems. 
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FIG. 1. Electron densities ph(r) and pc(f) for N interact- 
ing electrons in the disk geometry. The dashed lines show the 
Laughlin state for a short-range interaction, solid lines show 
the exact ground state for the Coulomb interaction. 



TABLE I. Normalized overlaps of the Laughlin ground 
state \Pl with the exact Coulomb ground state ^c, and the 
Coulomb ground state in the restricted Hilbert space for 
jV electrons. For the single-particle angular momentum 
m is restricted to m < m^ax, where mj^ax is for Here 
W is the weight of the components in the exact with 



a For N — 12 "unrestricted" overlap was obtained by an ex- 
trapolation, using the results for m^ax < "imax + 2 . 
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0.993606 


0.990997 


3.81 


4 


0.983108 


0.978804 


5.20 


5 


0.988168 


0.985047 


4.10 


6 


0.986102 


0.981767 


4.61 


7 


0.964758 


0.956392 


6.92 


8 


0.968652 


0.962187 


6.49 


9 


0.972992 


0.966526 


6.61 


10 


0.971415 


0.964413 


7.08 


11 


0.966118 


0.958138 


7.67 


12 


0.961240 


0.9528 a 


7.7" 
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For N > 10 the Laughlin state develops a constant 
density plateau pl ~ g in the interior of the disk. This 
is the precursor of the translationally and rotationally 
invariant v = ^ bulk FQH condensate. Surprisingly a 
second oscillation in the Coulomb pc {r) becomes evident 
for N = 10; the second oscillation shifts to larger r as 
N is increased, following the first oscillation. The sec- 
ond oscillation is clearly caused by the long-range part 
of the Coulomb interaction, in other words, by the Hal- 
dane pseudopotentials Vi with I > 3 (because of Fermi 
statistics, / must be odd for spin-polarized electrons). In 
view of large overlaps, one may want to describe the edge 
charge density oscillations in m terms of excitations 
of \&l by the long range part of the Coulomb interaction. 
Since the charge and the angular momentum M of the 
total system are fixed, a relevant elementary excitation to 
consider is a quasiparticle-quasihole bound exciton [po[ . 
The excitations that do not conserve m are not relevant 
because they alter M. A variational wave function that 
admixes one bound exciton into reproduces charge 
density oscillations in pc(r) for N = 10 to 12 reasonably 
well pl| . The deviation of the variational p(r) from the 
exact pc{ r ) increases, however, as N is increased from 10 
to 11 to 12, which indicates that more than one bound 
exciton is needed for larger systems. 

We have obtained a good phenomenological fit of the 
exact pc( r ) dependence using an analytical expression: 

p(r)=p [(erff+l)/2]{l 



where po = h for v 



/0ESP()Jo[<?ESp(y - »"o)]} ' ( 4 ) 
r is the distance 
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from the edge in units of do, erf (a;) is the Gauss error 
function, and Jo is the Bessel function of the first kind. 
The fitting parameters are the ESP oscillation amplitude 
at the edge p ES po, the ordering wave vector q ESPl and the 
shift of the first oscillation from the edge ?q. The fit 
is robust in the sense that the same fitting parameters 
also give good fits for N = 10 and 11, scaling r d g0 as 
(2mJ; ax ) 1 / 2 . A fit for N = 11,12 is shown in Fig. 2, 
where the values of the parameters are: p E spo = 0.519, 
wave vector q ESP — 1.570/£ , and f = 0.99. We be- 
lieve that the deviation of the fit from the exact density 
in the center of the disk for r < 1 is a finite-size effect. 
The parameters p E spo and f may be varied by 1% to give 
equally good overall fit. The wavelength A ESP = 4.002 can 
be compared to several natural lengths in the problem. 
It is significantly smaller than the Wigner crystal lattice 
constant a wc = (47rv / 3) 1 ^ 2 = 4.665, and, more relevant, 
is accurately equal to the distance between the crystal 
planes h wc = ( v / 3/2)a wc = (3irV3) 1/2 = 4.040. We also 
note that q ESF is close to the wave vector of the bulk 
magnetoroton minimum q MR p3). Similar length scales 
are provided by many other quantities, e.g.: the "ion di- 
ameter" 2i?o = v^o^r = 4.342, a little too large; square 
root of the "disk area per electron" (27rTO^ ax /A^) 1 / 2 = 
[67r(iV - Vj/N] 1 /' 2 = 4.157 for N = 12 is closer, but the 
Pc{t) data do not show such dependence on N, Fig. 2. 
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FIG. 2. Electron density pc{r) for N Coulomb-interacting 
electrons. The N = 11 density data (circles) is shifted radi- 
ally by r cdgc (7V = 12) - r cdgo (iV = 11) w 0.378. The dashed 
line is the fit of Eq. (4) with parameters given in the text. 

In the disk geometry rotational symmetry (correspond- 
ing to translational invariance along a linear edge) of the 
chiral basis orbitals |m) ensures that azimuthal density 
modulations in p(9) are not possible. A broken symmetry 
ground state is possible only as a superposition of several 
rotationally invariant degenerate ground states at differ- 
ent total M. Such ground state would be pinned by disor- 
der and QHE would not be observable in experiments. A 
charge density wave (CDW) order is possible in the radial 
direction, however; in the thermodynamic limit |23f| , the 
edge confinement breaks the translational invariance and 
"excites" radial CDW. Because of the long-range nature 
of the Coulomb interaction, the radial CDW is expected 
to propagate deep into the interior of the electron system 
and form a locally anisotropic edge striped phase (ESP). 
Since there is abundant experimental evidence that FQH 
edges are conducting, ESP is a conductor along the edge, 
and an insulator in the perpendicular direction, in the 
sense that in a QH system radial perturbation will induce 
only azimuthal Hall current. Thus, ESP can be thought 
of as a quantum smectic liquid crystal phase |24[ at the 
edge of the system. 

Extrapolating to the thermodynamic limit, and setting 
Pespo = | j 9esp = f j the following simple expression de- 
scribes well the electron density at the v = | FQH edge: 

p(y) = ~(erfy + 1){1 + ~ J [ |(y - 1)]} , (5) 

where y is the distance from the edge. In the disk geome- 
try the azimuthal transport current density Jg oc d r p(r). 
Differentiating Eq. 6 we obtain edge current density 



^(yjoc^e-^l + ljof^y-l)]} 

-^erfy+lJJ^fo-l)]. 



(6) 
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FIG. 3. Electron density and transport current density 
(arb. units) vs. distance from the v = i FQH edge. Density 
pcij) for 12 electrons is also shown. The analytical expres- 
sions Eqs. 5 and 6 are inferred from the numerical data. 

In practice, the first term with exp(— y 2 ) vanishes after 
the first oscillation. The phenomenological expressions 
Eqs. 5 and 6 are plotted in Fig. 3. Large y asymptotic 
behavior of both Bessel functions is cx y" 1 ^ 2 cos(7ty/2) 
(with a \ phase shift). Thus, should Eqs. 5 and 6 prove 
applicable to linear edges in large systems, the ampli- 
tudes of the charge density and the transport current 
oscillations both decay only as a power law of distance 
from the edge, specifically as y~ x ^ 2 ■ We note that the 
cx y~ x l 2 asymptotic behavior is not well established by 
our numerical results. 

There are likely to be several important implications of 
the observed slow-decaying charge density oscillations at 
FQH edges. Note first that they do not affect topological 
QHE properties, such as the values of the quantized a xy 
and the charge of the bulk Laughlin quasiparticlcs . 
On the other hand, dynamic properties, where edge exci- 
tations at finite frequency and wave vector are important, 
are likely to be affected. For example, since edge excita- 
tions are expected to be transverse charge density waves 
H , it seems that the presence of ESP at QH edges should 
affect spectrum of edge excitations Js) . The striped phase 
may be the reason for the reported experimental depen- 
dence of velocity of propagation of edge magnetoplas- 
mons on frequency pq |. When there are several edge 
modes present, a striped phase extending hundreds of 
£q from the edge is likely to couple various modes, and 
this may explain the failure to observe multiple and/or 
counterpropagating modes in experiments [^5| . 

A complementary field-theoretic approach for the dy- 
namics of an edge channel treats it as a Luttinger liq- 
uid, focusing on the asymptotic long-distance, low-energy 
physics ^,0. It has been argued that many expected 
power law behaviors in this limit are completely deter- 
mined by the quantized value of a xy , and would not be 
affected by the specifics of the edge structure. The map- 



ping of a FQH edge onto the Luttinger liquid is based 
on rather general considerations, provided there is no 
cx 1/r interparticle interaction. In this work we specifi- 
cally show that the long-ranged true Coulomb interaction 
leads to formation of a striped phase at the v — i FQH 
edge, calling into question the fundamental assumption 
that all important physics at the edge is one-dimensional. 
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